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PURPOSE 


The purpose of this document is to specify the equations 
necessary to perform the guidance, navigation and control onboard computa- 
tion functions for the space shuttle orbiter vehicle. This equations 
document will: (1) establish more specifically, than on a functional 

level, the GN&C computational requirements for computer sizing, (2) provide 
GN&G design equations specification to develop demonstration software for 
hardware feasibility testing, and (3) define the hardware interface 
requirements with the GN&C subsystem software. The document will provide 
a standard of communication of information concerning the GN&C equations, 
and will provide a means of coordination of GN&C equation development. 
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SCOPE 


This is Volume II of the document which defines the Guidance, 
Navigation and Control (GN&C) design equations sequencing and interfaces 
for the computations required in the GN&C Subsystem for all mission 
phases of the Shuttle Orbiter flight. The equations are intended to 
satisfy the functional requirements specified in Reference a. This 
document will describe in mathematical, logical, and operational language 
all the details necessary to initiate and carry out the design of the 
required computer modules (subprograms) for the GN&C functions. 

The document will be organized into six volumes. Volume I 
contains Sections 1 through 8, which provide introductory information for 
the document. Volume II contains the current detailed equations for the 
preflight, boost, separation, orbit insertion and ascent abort phases of 
the Orbiter operation. Detailed equations for orbital operations of the 
Orbiter, which include the orbital coast, orbital powered flight, rendez- 
vous, station keeping, docking and undocking, and docked operations phases, 
are presented in Volume III. Volume IV contains the current detailed 
equations for the deorbit and entry, transition, cruise and ferry, approach 
and landing, and horizontal takeoff phases of the Orbiter. Also, this 
volume will contain the equations for communications and pointing functions 
and the failure detection function. For this issue of the document, only 
Volumes I through IV are being published. In future issues. Volume V will 
contain the detailed flow diagrams for the equations. For the initial 
issue, the flow diagrams for the approved equations are included with the 
equations in Volumes II, III or IV. For future issues, the constants used 
in the equations will be summarized in Section 11 and the GN&C parameters 
and variables which can be entered or called via -the keyboard will be 
enumerated in Section 12. These two sections will be contained in 
Volume VI of the document. 
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APPLICABILITY 


This document is applicable to the Guidance, Navigation, and 
Control (GN&G) Subsystem of the Electronics System of the Space Shuttle 
Orbiter Vehicle. It is applicable to the definition of the shuttle 
computational requirements for the subsystem listed above. It is applicable 
to the Phase B and Phase C subsystem development. It defines the Manned 
Spacecraft Center Guidance and Control Division inhouse study baseline 
equations design. 
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9. DESCRIPTIONS OF EQUATIONS 

The detailed equations for the GN&C functions are defined in this 
section. The organization of this section is tentative and will be 
modified so as to present the equations as they are designed in as clear 
a fashion as possible. As an introduction to each major subsection 
(usually a mission phase) , the general GN&C software functions to be 
implemented will be identified and, where appropriate, a conceptual 
discussion and top level flow of the computations, inputs and outputs 
will be included in order to understand and summarize what is to be 
covered. This should be an order of magnitude less detailed than the 
flow diagrams of the equations which come later. 

A GN&C Equation Submittal sheet will introduce each of the GN&C 
equation submittals and summarize the GN&C functions, and identify the 
source and NASA contact for each. 

The detailed data to be presented for each GN&C function within each 
of the major subsections (usually a mission phase) is summarized below. 
Although items 6 through 10 are to be referenced only in the equations 
document, they are required submittals before the equations can be 
approved and finalized for flight software development. 

1. Functional Requirements 

The specific functional requirements (from the GN&C 
Software Functional Requirements Document) which are 
satisfied by the equations should be identified. 

2 . Functional Diagram 

A brief functional explanation and description of the 
overall concept and approach. A functional block 
diagram should be used where clarity is enhanced. 

Inputs, outputs, and Interfaces will be provided. 

3 . Equations and Flows 

Detailed equations and a descriptive text which guides 
the reader through the flows of Section 10 should be 
provided. The minimum frequency of the computations 
shall be specified and rationale given or referenced. 
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4 . Coordinate System 

The coordinate systems used shall be defined. 

5 . Cons tants /Variables Summary 

Constants and variables shall be summarized in tabular 
form with the following information: 

a. Variables /constants symbols and definitions 

b . Units 

c. Allot/ able quantization 

d. Range of values 

6. FORTRAN Coding 

The FORTRAN coding of the function for verification using 
the Space Shuttle Flight Simulation (SSFS) will be 
referenced . 

7 . Simulation 

The SSFS specifications, description and user's guide 
used to verify each GN&C function will be referenced. 

8. Testing 

Test plans and test results will be referenced. 

9 . Derivation 

The mathematical derivation of the equations including 
all mathematical assumptions shall be referenced. 

10. Assumptions 

The following will be referenced: 

a. Avionics baseline system assumed 

b. Reference missions assumed 

c. Vehicle mass properties assumed 

d. Propulsion models assumed 

e. Environment models assumed 

f. Error models assumed 

The major subsections of this section are identified and partially 
expanded in the following. 
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9.1 PKEFLIGHT 


The following GN&C software functions are envisioned for this mission 

phase: 

1) Generate a reference mission profile. Specifically, 
perform targeting for entire mission during the prelaunch 
phase to insure that mission objectives will be met. 

Target the launch phase to insure that launch abort 
targeting, rendezvous targeting, and deorbit targeting 
solutions are consistent with vehicle consumables, 
characteristics, etc., and ground related activities 

and locations . 

2) Schedule sensor calibrations (if required in flight)^ 
consistent with sensor status and mission plan. 

3) Schedule sensor alignments consistent with sensor status 
and mission plan . 

4) Calculate initial inertial state prior to launch. 

5) Transform launch guidance target parameters to proper 
time dependent form for launch. 

6) Navigate vehicle and target to launch time. 

7) Perform prelaunch checkout of all control hardware using 
autopilot and manual controller inputs. 

8) Perform automatic calibration of sensors (i.e., accelero- 
meters and gyros) or use predicted values based on time 
histories and compute compensation values prior to launch. 

9) Provide automatic inertial reference alignment prior to 
launch . 

9.1.1 GN&C System Initialization, Monitor, Test and Checkout (TBD) 

9.1.2 Prelaunch Targeting (TBD) 

9.1.3 Sensor Calibration and Alignment (TBD) 
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9.2 BOOST PEASE PROGRAM 


The Orbiter Boost Phase Program will monitor the Booster performance 
during the first-stage boost period. State vector and attitude data will 
be maintained and updated periodically so that the Orbiter will be ready 
for staging at any time. Figure 1 represents a functional flow of 
guidance, navigation and control computations during Boost. Redundancy 
internal to the orbiter is undefined at this time and is not shown. 
Navigation and Guidance computations will be performed throughout the 
Boost Phase in order to monitor the behavior of the vehicle, provide data 
to the abort decision process, and to prepare the orbiter for staging 
and orbit insertion. 

The guidance computations will include steering, throttle and engine 
cut-off time computations, although no takeover of Booster control is 
currently anticipated. The equations used for the guidance functions 
would be used only to detect malfunctions and would necessarily be the 
same equations used by the Booster. 

The Boost Phase Program would be terminated at staging and replaced 
by the Orbit Insertion Program or, following an abort decision, by the 
Ascent Abort Program. 
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Initialization (Prelaunch) 


9.2 BOOST PHASE PROGRAM (continued) 
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9.2 BOOST PHASE PROGRAM (continued) 


SPACE SHUTTLE 

GN&C SOFTWARE EQUATION SUBMITTAL 

Software Equation Section: Rapid 9 Real-Time State Advancement During 

Specific Force Sensing 

Submittal No. 19 

Function: Powered Flight Navigation 

Module No. ONI Function No. 1 (M SG 03690) 

Submitted by: W. M. Robertson Co. MIT No. 5-71 

Date : Feb 1971 

NASA Contact: J . Suddath Organization: GCD 

Approved by Panel III: |£.T. Ce£. Date: 3// 0 / 7 / 

f & 

Summary Description: Provides the capability to quickly and continuously 

update the spacecraft state vector in. real time, incorporating into the 
state update the real-time measurements from the accelerometers . 

Shuttle Configuration: Not provided explicitly. However. AV's from • 

inertiallv stable accelerometers are implied. 

Comments : 

(Design Status) 

(Verification Status) 

Panel Comments : 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing 

1. INTRODUCTION 


The routine for Rapid Real-Time State Advancement during 
Specific Force Sensing provides the capability to quickly and contin- 
uously update the spacecraft state vector in real-time (especially 
during periods of high computer activity), incorporating into the 
state update the real-time measurements which are being provided 
by the spacecraft’s accelerometers relative to the specific (non- 
free-fall) forces to which the spacecraft is Deing subjected, such as 
thrust and aerodynamic effects. 

The routine is merely a low order numerical integration of 
the differential equation satisfied by the spacecraft state vector in 
which only the central force field term and the largest earth oblate- 
ness term ( J 2 ) are included, but in which provision is made to 
incorporate the specific force data from the accelerometers. 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing ‘ (continued) 

NOMENCLATURE 


COS <j> 
J 2 

£<t) 

r 

r (t) 

r E 

v(t) 

At 


Av (At) 
-^neas 



1 

— z 


Cosine of colatitude of spacecraft position 
Earth Oblateness Constant 
Gravity vector at time t 
Magnitude of r 

Spacecraft position vector at time t 
Earth's mean equatorial radius 
Spacecraft velocity vector at time t 
Time interval 

Velocity vector increment measured by the integrating 
accelerometers during the time -interval At 

Earth's gravitational parameter 

Unit vector in the direction of r 

Unit vector in direction of earth's north polar axis 
expressed in reference coordinates 
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9.2.1 


Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 

2. FUNCTIONAL FLOW DIAGRAM 

The general functions performed by the routine are outlined 
in Fig. 1. 

The routine is called repetitively during specific force 
sensing periods once every At seconds, where the time-interval At 
is chosen so that the procedure will be compatible with the basic 
guidance and control computation cycles. 

Because the accelerometers are of the integrating type, the 
accelerometer data is output in the form of velocity increments 
( Av ) over the repetitive computation time interval At, The 
computations of this routine are therefore in terms of increments of 
velocity rather than instantaneous accelerations. The increments 
are initially sensed in Inertial Measurement Unit or Stable Member 
Coordinates but are transformed into reference coordinates before 
use in this routine. 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 


Enter from Guidance and Control 



Exit to Guidance and Control 


FUNCTIONAL FLOW DIAGRAM 











9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 

3. PROGRAM INPUT -OUTPUT 

The required inputs to the routine consist of the most recent 
state vector, the time interval through which the advancement of that 
state vector is desired, the velocity increment due to the specific 
forces which was measured by the accelerometers over this time 
interval, and the gravity vector at the most recent position vector. 
The outputs of the routine are the new ( i. e. advanced ) state vector 
and the gravity vector at the new position vector, which will be used 
as inputs to the next computation performed by the routine. 

Input Parameters 


System 


E 


k E 


1 

-2 


Earth's gravitational parameter 
Earth's mean equatorial radius 
Earth oblateness constant 

Unit vector in direction of earth's north polar axis 
expressed in reference coordinates 


Real-Time 


r(t), v(t ) 


Most recent spacecraft state vector ( at time 
t) 


At : Time interval through which the most recent 

spacecraft state vector is to be advanced 

A v (At) : Velocity vector increment which was measured 

— m “3.S 

by the integrating accelerometers during the 
time interval from t to t + At 
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g(t> 


Gravity vector at the most recent position 
vector r(t). 


Output Parameters 


New (advanced) spacecraft state vector 
( at time t + At ) 

g ( t + At ) : Gravity vector at the new position vector 

r (t + At). 


r (t + At ), 
v ( t + At ) 
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9-2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 


4. DESCRIPTION OF EQUATIONS 

The position and velocity vectors are advanced through a 
time intervalAt by the formulae: 


r(t At) = r ( t) - 1 - At J~v(t) x ^-|[(t) x - ^^ meas .. 


] 


At 


v(t x At) = v(t) + yi| £(t x At) + g(t) 


£ 


+ Av (At ) 
-meas 


where Av mgag (At) represents the velocity vector increment mea- 
sured by the integrating accelerometers during the time interval At. 
The formulae represent a second -order technique of numerically in- 
tergrating the differential equations of motion 

j2 

- d - r(t) = £ 
dt 2 

where g is a simplified form of the gravity vector at the position r. 
The technique requires only one evaluation of £ for each time-step 
At, which is of course a definite advantage in reducing overall compu- 
tation time. 

The gravity vector £ at the position r is calculated from the 
expression: 



1 

2 

“ — 

1 

M E J ^3 

r E 


(1-5 cos 2 d>) 1 x 2 cos <4 1 

l 

£ = ‘ ~2 S Ar 2 J 2 

r 


— r — z 

( 

r l 

i 


_ rr 

J 


where 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 

1 = unit vector in the direction of r 

1 ^ = unit vector in the direction of the earth’s polar 

axis expressed in reference coordinates 

cos d — 1 • 1 

— Z 

and 

r^- = the earth's mean equatorial radius. 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 


5. DETAILED FLOW DIAGRAM 


This section contains a detailed flow diagram of the routine 
for Rapid Real-Time State Advancement during Specific Force Sens- 
ing, 
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9.2.1 Rapid Real-Time State Advancement During 
Spacific Force Sensing (continued) 


Enter from Guidance and Control 



Exit to Guidance and Control 


Figure 2 DETAILED FLOW DLAGEAM 
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9.2.1 Rapid Real-Time State Advancement During 
Specific Force Sensing (continued) 


6 • S UPPLEMENTARY INFORMATION 

The formulation was originally developed for the Apollo 
Guidance Computer by Bairnsfather (1963 ), who also gave an error 
analysis of the computation when performed in eight, nine, and ten 
decimal digit precision ( with addition instructions both truncated and 
rounded). It should be noted that Bairnsfather 's analysis shows that 
a time -step of two seconds yields the most accurate results in terms 
of both position error magnitude and velocity error magnitude for 
eight decimal digit precision with rounding over a thrity-five minute 
integration interval (100 ft, and 0,2 ft/ sec, respectively). 


References 


1. Bairnsfather, R., 1963, Navigation Algorithms for Non- 

Freefall Flight , MIT/IL Space Guidance Analysis Memo 
#47, June 1963, 


The Apollo Guidance Computer has somewhat over eight decimal 
digit precision. The guidance programs written for it have had 
(for a variety of reasons) a guidance cycle period of two seconds. 


9 . 2-14 


9.2.2 Boost Monitor (TBD) 




9.3 SEPARATION 


The Separation Phase covers the portion of the orbital mission 
between the preparation for staging (during Boost) and the ignition 
of the Orbiter main engines. The maneuvers required during separation 
will be highly dependent on the mechanics of separation and the restric- 
tions on relative vehicle motion due to impingement problems. Control 
functions during separation are greatly dependent on such open items as 
the vehicle configurations before and after separation, the separation 
mechanism employed, the plume impingement criteria, and the coast time 
limits . 

The general GN&C software functions to be performed during the 
Separation are the following : 

1. Continue powered flight navigation. 

2. Command Staging sequence. 

3. Compute and command ECS control during separation. 

4. Perform attitude control during coast. 

The following provides a preliminary description of the Shuttle 
dynamics during nominal and abort staging maneuvers. 

"Preliminary Shuttle Dynamics Analysis," MSC 
Internal Note No. 70-PM-143 , 23 September 1970. 

9.3.1 Attitude Control (TBD) 
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9.4 ORBIT INSERTION PROGRAM 


The Orbit Insertion Phase extends from separation from the Booster 
through the second-stage boost until the desired orbital velocity has 
been achieved or an abort decision has been made. Presumably, data for 
initialization has been computed previously during the Boost and 
Prelaunch Phase. The Orbit Insertion Program is terminated by the start 
of the Orbital Coast Program or the Ascent Abort Program. 

9.4.1 Navigation (Same as Rapid Real-Time State Advancement 
During Specific Force Sensing) 
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SPACE SHUTTLE 


GN&C SOFTWARE EQUATION SUBMITTAL 


Software Equation Section: Powered Ascent Guidance Submittal No. 4^ 

Function: Provide Guidance Commands During Orbit Insertion 

Module No. 0G1 Function No. 1 (MSC 03690) 


Submitted by: J. P. Higgins Co. MIT No, 1-71 

Date : Jan 1971 

NASA Contact: W. Peters Organization: GCD 

Approved by Panel III: Date: Zl/ ZV/ 7 / 

Summary Description: The orbit insertion maneuver is controlled so as to 

achieve a specific orbit. To control the ascent to the position and veloc ity 
.constraints, an explicit guidance law employing linear control is used. 

.The radial and cross-range position and velocity constraints are controlled 
by their acceleration components. Downrange is controlled b y time of 
thrust termination. ~ 


Shuttle Configuration: Equations are presented symbolically so that the 

configuration need not be specified. 

Comments : 

(Design Status) 

(Verification Status) 

Panel Comments : 
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9.4.2 Powered Ascent Guidance 
1 INTRODUCTION 

The objective of the Ascent Guidance Program is to provide 
guidance during the orbit insertion maneuver. This maneuver is con- 
trolled so as to achieve a specified orbit, which can be defined in 
terms of a desired velocity vector, a desired radius and a desired 
orbital plane at injection. To control the ascent maneuver to these 
position and velocity constraints, an explicit guidance law employing 
linear control is used. The radial and cross-range position and 
velocity constraints are explicitly controlled by two acceleration 
components, ( a_ r ) and ( a^-^} 5 that vary linearly with time in the 
radial and cross-range directions. The third component of accelera- 
tion is determinate. The horizontal, down-range velocity constraint 
is controlled by terminating thrust at the proper time. 
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9.4.2 Powered Ascent Guidance (continued) 


NOMENCLATURE 


A A A A 
11 ' 12 ' 21 ' 22 


TL 


*Tr 

a Ty 

a Tz 

B 

C 

D 

^eff 


i-d 

^-q 


— r 


l 

— z 


r 

r 

r 

r 


linear -control coefficient 

parameters needed to compute linear -control 
coefficients 

estimated thrust acceleration 
thrust-acceleration limit 
required radial acceleration 
required cross -range acceleration 
required down-range acceleration 

linear -control coefficients 

effective gravity 

unit vector in the required thrust direction 

unit vector normal to the desired orbital 
plane 

unit local-vertical vector 
unit cross -range vector 
unit down-range vector 
vehicle -position vector 
radius 

radial-velocity magnitude 
desired radius 

desired radial -velocity magnitude 

current time associated with r and v 

time at which the thrust direction command 
is issued 
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9.4.2 Powered Ascent Guidance (continued) 


At 


t. 


g° 


V 


V , 


V 


V 

-g 


y 

z 

*d 

?d 


X 

u 

y. 


guidance -cycle time step 
ignition time 

time -to -go until thrust termination 

time -to -go until acceleration limit is reached 

vehicle velocity vector 

desired velocity magnitude 

effective exhaust velocity 

velocity -to -be -gained vector 

cross -range distance measured from the 
desired orbital plane 

cross -range velocity magnitude 

down-range velocity magnitude 

desired down-range velocity magnitude 

desired flight path angle 

desired inclination angle 

throttle setting - fraction of full throttle 

earth's gravitational constant 

desired longitude of the descending node 
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9.4.2 Powered Ascent Guidance (continued) 


2 FUNCTIONAL FLOW DIAGRAM 


The chronological sequencing of functions performed by the 
Ascent Guidance Program is described in this section and illustrated 
by the functional flow diagram in Figure 1 . 

The Ascent Guidance Program is entered directly from the 
Boost Monitor Program when separation of the orbiter stage from 
the booster stage is sensed. The first function performed is to verify 
ignition time and to adjust it, if needed. Next, the desired target 
parameters are initialized and the initial thrust direction is computed. 
After these initialization functions are completed, the program cycles 
through the powered flight navigation routine until ignition occurs. 
Once ignition occurs, steering is initiated and the program cycles 
through the navigation routine and the steering routine every time 
step until thrusting is terminated. 
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9 . 4.2 


Powered Ascent Guidance (continued) 


ENTER 



EXIT 


figure 1 FUNCTIONAL FLOW DIAGRAM 
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9.4.2 Powered Ascent Guidance (continued) 


3 PROGRAM INPUT -OUTPUT 

The required inputs to the Ascent Guidance Program are the 
ignition time, an estimate of time-to-go, and a set of parameters de- 
fining the injection velocity vector, radius magnitude, and orbital 
plane. The output of this program is the thrust direction required at 
the present time to achieve the desired orbit and the throttle setting 
required to avoid exceeding the acceleration limit. 

Input Parameters 


go 

t. 

ig 


v d 

'Id 

l d 


desired radius 

estimate of the time -to -go until thrust termination 
ignition time 

desired velocity magnitude 

desired flight path angle 

desired inclination angle 

desired longitude of the descending node 


Output Parameters 
i^ required thrust direction 

A required throttle setting 
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9.4.2 Powered Ascent Guidance (continued) 


4 DESCRIPTION OF EQUATIONS 


The Ascent Guidance Program has two major phases, the 
pre -ignition phase and the powered ascent phase. Section 4. 1 describes 
the equations required for the pre-ignition phase and Section 4. 2 de- 
scribes the equations required for the powered ascent phase. The 
computational sequences are illustrated in Figs. 2 and 3. 

4. 1 Pre -Ignition Phase 

Figure 2 outlines the computational sequence of the pre-ignition 
phase. After ignition time (t. ) has been verified, the set of para- 
meters defining the desired injection velocity vector, radius magnitude, 
and orbital plane are input. Then the following thrust parameters are 
initialized: 


TL 

v 


= ( tbd ) 

= ( TBD ) 
= ( TBD ) 
= ( TBD ) 


Initial thrust filter parameters = ( TBD ) 


Next, the time -to -go until thrust termination ( t ) is set to the nom- 
inal value and the desired injection target parameters are initialized 
as follows. 


go 


= (TBD) 


= v ri sin 


V, cos (v ,) 
a ■ a 


i 

-q 


-sin U d ) sin ( 
sin (l d ) cos ( j 5 d ) 
cos ( l d ) 
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9.4.2 Powered Ascent Guidance (continued) 


Fin all y, the initial required thrust direction is computed by cycling 
to point (A) of Fig. 3. After these pre-ignition computations are com- 
pleted, the program cycles through the navigation subroutine until 
ignition is sensed. At this time, the program enters the powered as- 
cent phase. 

4. 2 Powered Ascent Phase 

The computational sequence during the powered ascent phase 
is outlined in functional form in F ig. 3 and the related equations are 
described below. These equations are recomputed at each time step 
during the powered maneuver unless otherwise stated. 

Thrust Filter and Throttle Control ( TBD ) 

The purpose of the thrust filter computations is twofold. 

First, the incremental velocity changes measured by the accelero- 
meters include noise. As a result, the thrust filter is necessary to 
accurately estimate the present thrust acceleration ( a^ ) and the mass 
to mass-flow-rate ratio (r). Secondly, since thrust acceleration is 
limited, it is necessary to predict the thrust acceleration over the 
next time step in order to determine if the engines need to be throttled 
down so as not to exceed this limit. The output of the thrust compu- 
tations will be the following. 

( 1 ) a j = estimated thrust acceleration 

(2) r = mass to mass -flow- rate ratio 

( 3 ) ^ = throttle setting 


Local Position and Velocity Computations 

After the thrust filter computations are completed the follow- 
ing parameters are computed : 
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9.4.2 Powered Ascent Guidance (continued) 


i = unit ( r) 


i = unit (i„xi ) 
— z — r ~q 


i = i x 1 
-y -z -r 

r = 1 r I 


y = r sin 1 ( i • i ) 
J — r -q 


r = v • i_ r 


y = v • i 

_ - -y 

z = v • i 
— — z 


Effective Gravity 


Next, the effective gravity is computed for use in the velocity - 
to -be -gained computation, 

g eff = ~ u ‘ r ll x l\ I r 


Velocity -to -be -Gained 


The most recent estimate of the time-to-go is decremented 
by the guidance -cycle time step. 


t 


go 



At 


The velocity -to -be -gained is then computed by differencing the present 
velocity components with the desired velocity components and correct- 
ing for the average effect of gravity over the time-to-go. 


Vg = (r d - r) i_ r -y iy Mz d - z ) i_ z 
1 . 

V = V £ t 1 

-g -g 2 ^ eff go -r 
v = | v } 

g l-gl 

Time -To -Go Prediction 


The usual exponential equation 
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9.4.2 Powered Ascent Guidance (continued) 


go 


= t ( 1 - e 


-v /v 

rtf c. 


valid for a fixed-thrust rocket must be modified to compensate for 
g -limiting. First a t computation which neglects g -limiting is 
computed from a truncated series expansion of the previous equation 


go 


V , V 

= t _£ 


V 


2 v e 


Then the time-to-go until the acceleration limit is reached ( tj- ) is 
computed from the following equation 


t 


L 


= T - v 


/ a TL 


If tg Q is greater than t^, then t^ Q is recomputed to account for g- 
limiting by the following equation 

*80 = *L [ V g " V e log ( 1 - *l7 T) ] ' a TL 
Guidance Parameters 


The Ascent Guidance Program uses a linear control law which 
can be expressed by the following equations. 

a Tr = a T " B(t ' V " g e ff 


a Ty _ a T j_ C ^ D ft t Q ) J 


^ are the required radial and cross-range accelerations 
respectively, is the filtered thrust acceleration, -and (t - t ) is 


a_ and a 
i r 


the difference between the time when the linear coefficients are 
computed and the time when the thrust direction command is issued. 
When t becomes less than At seconds, the linear coefficients 

go 

are not recomputed; instead, the previous values are used to compute 
a— ^ and a__. The equations required to compute the linear coeffi- 
cients are as follows: 
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9.4.2 Powered Ascent Guidance (continued) 


A =[- A l2 (r d ' r ' ^ 'go >^22^ 

B= [ A 11 (r d - r ' * 'go > - A 21 (i- d 

c=[ a 12 (y+yt go ) - A 22 y] J 
D =[ - A u (y + yt go )+A 21 y J J 

where 

J = a 22 - a 12 a 21 ) 

The parameters A^. A l2 , A 21 , and A 22 are the coefficients of the 
simultaneous equations used to solve for A, B, C, andD. These para- 
meters are computed in one of three ways, depending on the current 

values of t and t . 

go L 

1 ) t „ < t T (acceleration limit will not be reached) 
go L 

A 11 = - v e lo S (1 ~ 'go lT) 

A 12 = - v e Tlo § (1 - 'go / T) - v e 'go 

A 2 1 V e (T -'go ‘'go /T)Xv e 'go 

A 22 = v e T(T ' 'go )l°g<l-t go It) 

j-v t ( T - i t ) 
e go 2 go 
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9.4.2 Powered Ascent Guidance (continued) 


2 ) 0 < t T < t ( acceleration limit will be reached ) 

L go 


A H = ' v e logd -t L /r) + a TL ( t gQ - t L ) 


A 12 = “ T log ( 1 - t T /t ) - v, t T + 


e ‘l 2 TL ' go 


( t -t T a ) 


A 21= v e {T - t go )iog(l-t L /r) + v e t L 

+ 2" a TL ( t go " *L )2 

A 22 = V e T(T ' t go )logd-t L /T)+v o t T (r + i-t T -t_ ) 


d ‘'T ' • o" U T L „. 

e Lu 2 L go 


+ - a TT (t 3 - 3 t t T 2 + 2t T 3 ) 

g TL go go L L ' 


3 ) < 0 ( acceleration limit has been reached ) 


A 11 " a TL t go 

A -I a t 2 
a 12 2 a TL go 

a 1 t 2 

21 2 a TL go 

A = — a t 3 
A 22 6 a TL go 


Required Thrust Direction 

The final computation performed by the Ascent Guidance 
Program is the required thrust direction. The radial and cross-range 
accelerations required to satisfy the linear control law are combined 
as follows: 

-Try Tr -r iy -y 


1 Try ! -Try | 
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9.4.2 Powered Ascent Guidance (continued) 


a Try is then compared with a T - If a T f s i ess than a Try’ the ac_ 
celeration vector a rp r y is reduced by the ratio a^/a^^. 

-Try = -Try a T / a Try 

If a T is greater than a,^ 4 the additional available acceleration is 

applied in the down-range direction. 

_ , 2 2 , 1 / 2 . ’ „ , 

a Tz ~ a T ' a Try ^ / sl S n ( z d " z 1 

The total thrust acceleration vector then becomes 

-T = -Try + a Tz iz 
and the required unit thrust vector is 

J d = unit (a T ) 
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9.4.2 Powered Ascent Guidance (continued) 

















9.4.2 Powered Ascent Guidance (continued) 


5 DETAILED FLOW DIAGRAMS 

This section contains detailed flow diagrams of the Ascent 
Guidance Program. 
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9.4.2 Powared Ascent Guidance (continued) 


ENTER 



CYCLE TO POINT 



of Fig. 5a 


Figure 4 DETAILED FLOW DIAGRAM, PEE -IGNITION PHASE 
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9.4.2 Powered Ascent Guidance (continued) 


Nominal Entry from 
Navigation (Every Guidance Cycle) 


Initial Entry from 
Pre -ignition Phase 



Figure 5a DETAILED FLOW DIAGRAM, POWERED ASCENT PHASE 
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9.4,2 Powered Ascent Guidance (continued) 



Figure 5b DETAILED FLOW DIAGRAM, POWERED ASCENT PHASE 
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9.4.2 Powered Ascent Guidance (continued) 



V 

Figure 5c DETAILED FLOW DIAGRAM, POWERED ASCENT PHASE 
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EXIT TO NAVIGATION 


Figure 5d DETAILED FLOW DIAGRAM., POWERED ASCENT PHASE 
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9.4.2 Powered Ascent Guidance (continued) 


6 SUPPLEMENTARY INFORMATION 

The guidance equations presented here were derived from 

explicit guidance laws developed by G. Cherry 1 . They were previously 

2 

used in the Apollo Lunar Ascent Guidance Program . Since the space 
shuttle acceleration is limited, however, a major modification of the 
equations was required. The derivation of this modification can be 
found in Appendix A of Reference 3. 

These equations have been used satisfactorily in the generation 
of various launch trajectories required for IMU performance analysis. 
In addition, a comparison was made with an optimal trajectory in one 
case. The results of this case showed the ascent-guidance equations 
to be near fuel optimal. Details of this comparison can be found in 
Reference 4. 

Future work should include development of a thrust estima- 
tion filter, throttling logic, and thrust tailoff compensation. Since 
they are dependent on detailed characteristics of the ascent engines, 
they have been omitted from this document. In addition, a compari- 
son of this scheme with other candidates, in particular the Saturn 
IGM equations, should be made. 
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Powered Ascent Guidance (con't) 
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9.4.3 Attitude Control (TBD) 
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9.5 ASCENT ABORTS 


Neither equations nor functions have been defined for Aborts during 
ascent. It is assumed here that the Abort Program will employ routines 
from other programs for most GN&C activities. The targeting would be 
performed by the Abort Program until the Orbiter has achieved a flight 
condition suitable for normal program sequencing for the return to Earth. 
The following functions are postulated for the Abort Program: 

1. Target the Orbiter to a parking orbit, a suitable entry 
condition, or an acceptable cruise condition. 

2 • Sequence other GN&C programs to perform maneuvers 
required during abort operations. 

9.5.1 Abort Decision and Mode Determination (TBD) 

9.5.2 Abort Targeting (TBD) 

9.5*3 Abort Guidance (TBD) 

9.5.4 Navigation (TBD) 

9.5.5 Attitude Control (TBD) 
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